Mon. Not. R. Astron. Soc. OOO,[TH5](2012) Printed 2 March 2013 (MN WI^ style file v2.2) 



The EoR Sensitivity of the Murchison Widefield Array 
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ABSTRACT 

Using the final 128 antenna locations of the Murchison Widefield Array (MWA), we 
calculate its sensitivity to the Epoch of Reionization (EoR) power spectrum of red- 
shifted 21 cm emission for a fiducial model and provide the tools to calculate the 
sensitivity for any model. Our calculation takes into account synthesis rotation, chro- 
matic and asymmetrical baseline effects, and excludes modes that will be contaminated 
by foreground subtraction. For the fiducial model, the MWA will be capable of a lAa 
detection of the EoR signal with one full season of observation on two fields (900 and 
700 hours). 
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1 INTRODUCTION 

The Cosmic Dark Ages and the Epoch of Reionization 
(EoR) remain a largely unexamined chapter of the history 
and evolution of the Universe. Observation of redshifted 
21 cm emission show s pr omise of probing the Eo R (see 
iFurlanetto et al" I l2006l and [Morales fc Wvithd [20T0I for re- 
cent reviews). Indeed, studies of the EoR and dark energy 
were rated with high priority by the 2010 Astronomy and As- 
trophysics Decadal Survey. Several ground-based radio ex- 
periments are currently under construction to probe the EoR 
through 21 cm power spectrum measurements, including 
LOFAR (LOw Frequency ArrajQ), PAPER (Precision Ar- 
ray for Probing the Epoch of ReionizatiorQ), and the MWA 
(Murchison Widefield ArrajO). 

The MWA is being built in the radio quiet Murchison 
Radio Observatory in Western Australia, and aims to mea- 
sure the EoR power spectrum via the 21 cm signal over a 
large range of redshifts. The originally planned MWA was to 
consist of 512 antennas, distributed over a circular region of 
radius 1.5 km (jLonsdale et al.|[2009l ). With current funding 
the instrument has been re-scoped to 128 antennas, but will 
have similar layout characteristics to the originally planned 
512 antenna array. A full description of the 128 antenna 
instrument is presented in lTingav et al.l (|2012l ). and a thor- 
ough description of the science capabilities will be presented 
in [Bowman et al. (2013). 

Here we calculate the MWA's expected sensitivity to 
the EoR signal using the physical antenna locations. A given 
baseline (the separation vector between any two antennas) 
is sensitive to a particular angular Fourier mode on the sky, 
so the baseline distrib ution is directl y related to the EoR 
sensitivity of an array ([Morales! l2005h . The MWA baseline 
distribution will have a dense core for EoR sensitivity and 
a smooth extended radial pr ofile for calibration and fore- 
ground subtraction purposes (|Bowman et al.ir2006l ). The lo- 
cations of the 128 antennas for the MWA were optim ized 
using the algorithm presented in lBeardslev et al] (|2012| ) and 
shown in Fig. [T] A table of the locations of all 128 antennas 
are available in the electronic supplement. 

We use a fiducial model to calculate the MWA's sensi- 
tivity and in attached tables provide the information needed 
to quickly apply any model. The EoR observing plan for the 
MWA is to track fields when they are above 45 degrees ele- 
vation and the sun and galactic center are below the horizon. 
Over an annual cycle, this yields a full observational season 
of 900 hours integration on a primary field and 700 hours 
on a second field. For the fiducial model, we find that with 
a full season of observation the MWA will be capable of a 
14cr power spectrum detection, along with constraints on the 
slope. 

Throughout this paper we use a ACDM cosmology with 
r2„ = 0.73, f^A = .27, and h = 0.7, cons istent with WMAP 
seven year results (jKomatsu et al.ll201ll ). All distances and 
wavenumbers are in comoving coordinates. 
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Figure 1. The antenna locations for the 128 antenna MWA. Po- 
sitions are measured relative to -26° 42' 4.396" Latitude, 116° 
40' 13.646" Longitude. The blue squares show the core 112 an- 
tennas which will be integrated for an EoR measurement. The 
solid red squares represent the outlier antennas used for solar 
measurements, but are not used for EoR measurements. While 
the antennas are indeed square, the squares shown here are not 
to scale. 



2 EOR SENSITIVITY 

The power spectrum measurement of the sky temperature 
is done in three dimensions (two angular directions, and the 
line-of-sight direction achieved through redshift), so we must 
find the uncertainty in each three dimensional voxel in cos- 
mological wavenumber (k) space, then perform a weighted 
average in spheric al bins to arrive at one-dimensional sen- 
sitiv ity (followi ng Moraled l2005l . iMcQuinn et al.ll2006l . and 
Morales fc Wvithe 201Q ). 

The fundamental visibility measurement of an interfer- 
ometer is done in (u, v, f) space, where u and v are the 
baseline coordinates (measured in wavelengths), and / is 
the frequency of the observation. The thermal uncertainty 
on the visibility measurement is given by 



(1) 



^ http: / /www. lofar.org/ 
^ http://eor.berkeley.edu/ 

http://www.mwatelescope.org/ 



where T^ys is the system temperature, is the effective col- 
lecting area per antenna. A/ is the frequency channel width, 
and r is the total i ntegration time for the mo de including re- 
dundant baselines (|Morales fc Wvithell20ld ). Observational 
parameters for our calculation are shown in Table[T] The sys- 
tem temperature is dominated by galactic fo reground emis- 
sion, and redshift dependence is discussed in lBowman et all 
(2006). Here we assume a constant system temperature over 
the observational bandwidth. 

To determine r, the integration time per (u, v, f) voxel, 
we use the surveyed antenna locations, and perform an aper- 
ture rotation for 3 hours on either side of zenith. We approx- 
imate chromatic effects by calculating the baseline migra- 
tion along the frequency dimension, then averaging. This 
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Table 1. Observational parameters for sensitivity estimation 



Parameters 



Values 



No. of antennas 
Central frequency 
Field of view 

Effective area per antenna 
Total bandwidth 

Tsys 

Channel width 
Latitude 

Primary Field RA 
Secondary Field RA 



112* 
158 MHz (z r. 
31° 
14.5 m^ 
8 MHz 
440 K 
40 kHz 
-26.701° 
6h 
Oh 



*Sixteen of the 128 antennas are not integrated for 
EoR measurements and are not included here. 
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Figure 2. Effective integration time per day per (n, v) cell includ- 
ing rotation synthesis for 128 antenna MWA at 2: 8 (A = 1.89 
m). The color scale units are the logarithm of effective seconds 
observed per day per angular mode, assuming six hours of inte- 
gration per day on one EoR field. Note that the total number of 
seconds observed per day is 21600, but an angular mode can be 
effectively observed longer due to redundant baselines. The most 
observed mode in this array is ~ 4 X 10^ sec/day. The uv cell 
size is dic tated by the size of th e instrumental window function. 
Following Bo wman et al. we used a cell size of (8.3 m)^. 

Data for this figure will be available in a machine readable table 
to easily plug into Equation [S] 



includes chromatic effects while avoiding a full covariance 
calculation (Hazelton, et al. 2012, in preparation). 

The sampling matrix for one day of observation on a 
single EoR field is shown in Fig. [2l The MWA will have a 
very dense, highly redundant uv core, with a smooth radial 
profile extending to 1.5 km. The large number of baselines 
in the core will beat down the thermal variance for those 
modes because the effective observing time is the sum of all 
the baselines observing the mode. 

The (u, v) coordinates map directly to the transverse 
cosmological wavenumber by the relation k± — 2nu/D, 



where D is the comoving distance to the observation. The 
observing frequency dimension maps to the line-of-sight di- 
rection, and must be Fourier transformed to the fc|| dimen- 
sion. Once these conversions are done, our data is in three 
dimensional fc-space, and we square to reach the power spec- 
trum. Propagating errors, the thermal uncertainty per k- 
space bin is given by 



C^(k) 



(2) 



The second term can be thought of as converting the uv 
bin size (As) to cosmological wavenumber space and has 
units of Mpc^. The third term converts the width of the 
observation from bandwidth to line of sight spatial extent 
and has units of Mpc s (for flat space the line-of-sight and 
transverse distances are equivalent), and t is the integration 
time for the k-space bin (in seconds). Inserting the values 
from Table [T] for all terms except the integration time gives 



C^(k) = ^-^^ ^ mK^Mpc^ 



(3) 



There is also a sample variance contribution to the un- 
certainty. Assuming the distribution is Gaussian, the sample 
variance per three dimensional voxel is given by the power 
spectrum itself (McQuinn et al. 2006). Combining the ther- 
mal and sample uncertainties gives the total variance per 
3D fc-space voxel 



<7|>(k)= P2i(k) + C^(k) 



(4) 



Because of the sample variance term, the calculated sen- 
sitivity of an array depends on one's choice of theoretical 
EoR model. While surveying the landscape of EoR models 
is beyond the scope of this paper, we have included a ta- 
ble of the effective seconds observed per day per {u, v) cell 
in the electronic supplement (data for Figure [2]). The sec- 
onds per day can be combined with the observing strategy 
to calculate the integration time per cell, r in Equation O 
and combined with the theoretical model in Equation |4] to 
accurately determine the sensitivity of the MWA for any 
proposed model. The coefficient values in Equation and 
the coordinates of the supplemental table can be scaled to 
different redshifts with ~ 5% error on the resulting sensi- 
tivity, or the antenna locations and synthesis rotation can 
be used to recalculate the integration time per bin using 
the supplemental table as a cross-check. In the remainder of 
this paper we use the fi ducial model of a fully neutral IGM 
iFurlanetto et al1l2006l 1^ as an example of how to accurately 
calculate the EoR power spectrum sensitivity. 

The underlying EoR fiuctuations are assumed to be 
isotropic, however velocity distortions will amplify the sig- 
nal in the line-of-sight direction on relevant large scales 
(jBarkana fc Loebll2005l l. For our fiducial model this angular 
dependence is given by P2i(k) = (l-|-2/i^ +/i'*)P2i(fc), where 
/i = fc||/|k|. This effect depends on whether dark matter or 
ionizing sources are sourcing the fiuctuations. Throughout 
reionization both sources will be relevant and the above ex- 
pression will depend on the cross-power spectrum between 
the fiuctuations. Because our fiducial model is a fully neu- 
tral IGM we can use this simplified relation. In addition the 



Available online at www.astro.ucla.edu/ sfurlane/21cm_pk.htm 
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Figure 3. Estimated power spectrum sensitivity to EoR 
signal per voxel for the MWA. The quantity plotted is 
log^Q(P2i (k)/(a-p (k)) for a two dimensional slice of the three di- 
mensional data cube with 900 hours of integration. The white 
curved lines show the bin edges used for the one dimensional plot 
(Fig. |4|| . The data below the horizontal dashed line and to the 
right of the diagonal dashed line will be contaminated by fore- 
grounds. Only data within the EoR window (the upper left) is 
used to calculate the sensitivity in Fig. (4] For reference, the cor- 
responding baseline lengths are given on the top axis. 
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Figure 4. Estimated ID sensitivity for the MWA, for various in- 
tegration scenarios. The dotted bl ue line is the theoretica l spher- 
ically averaged power spectrum llFurlanetto et aLll2006h . where 
A2j(fc) = P2i(fc)fc3/(2jr2r2) and To = 28[(1 -|- ^)/10]l/2 mK 
26.6 mK. The several step functions represent the uncertainty 
per bin, with the edges of the steps corresponding to the edges 
of the bins when averaging (white curves in Fig. I3J. Single field 
observations are shown for 450 hours (dash-dot green) and 900 
hours (dashed red) of integration. The solid black line corresponds 
to 900 hours on a primary EoR field, combined with 700 hours on 
a secondary field. This averaging excluded any data that would 
be contaminated by foreground subtraction (below the horizontal 
line, and to the right of the wedge in Figure |3j. 



MWA is sensitive to much smaller fcx modes compared to 
fell modes, so this effect is significant for the dark matter 
sourced fiducial model. 

Figure [3] shows the signal to noise per voxel in a slice 
through the 3D fc-space. At low k, a large signal and dense 
baseline distribution result in a signal to noise approaching 
1. Moving up in fc||, the signal diminishes, but the baseline 
density remains constant, so the sensitivity drops relatively 
slowly. Moving up in k±, however, both the signal and the 
baseline density drop, resulting in a more drastic drop in 
sensitivity. 

Foreground contamination limits the observability of 
the EoR. Fortunately, the contamination is localized in 3D 
fc-space, leaving a relativel y uncontaminate d EoR window 
l|Vedantham et all 2012: Mo rales et alll2012l '). The spectrally 
smooth foregrounds are fit to low order polynomials over 
the full 30.72 MHz instrumen t bandwidth, contam inating 
low line-of-sight wavenumbers (|Bowman et al.l [20091 ). How- 
ever, an individual observation is limited to ~ 8 MHz due 
to cosmic evolution, so only our A;|| = bin will be con- 
taminated. This exclusion zone is shown in Fig. [3] by the re- 
gion below the horizontal white line. In addition, mode mix- 
ing effects will throw pow er higher in fc||, cr e ating a wedge 
shap e of cont amination (|Datta et al.l l201ll : [Morales et al.l 
|2012| : [Xrott et al. 2012). The location of this contamination 
is indicated in Fig. [3] by the region below the diagonal line. 
The 'EoR window' is to the left of the diagonal line and 
above the horizontal line. In this calculation we only use 
modes within the EoR window. 

The next step is to perform a weighted average to con- 
dense the three dimensional data into a one dimensional 



power spectrum. The underlying power spectrum is ex- 
pected to be isotropic, so averaging in spherical shells of 
constant |k| is appropriate. As discussed earlier, the velocity 
distortion terms cause the power spectrum to be anisotropic, 
but can be remedied by dividing the signal and noise by the 
angular dependence, (1 -I- 2/x^ -I- /x*) in our case. Then voxels 
within a constant k shell have the same power spectrum sig- 
nal, and can be averaged weighting by the uncertainty per 
voxel. 

Figure |4] shows the sensitivity of the MWA to this EoR 
power spectrum. The theoretical one dimensional spherically 
averaged power spectrum (dotted blue line) and the uncer- 
tainty per bin (various step lines) are plotted. The uncer- 
tainty is plotted as a step function to show the binning used 
in the spherical average with the edges of the steps corre- 
sponding to the white curved lines in Figure |3l 

The uncertainty is shown for 450 and 900 hours on one 
field, as well as a two field observation with 900 hours on 
one field (RA = e*") and 700 hours on a second (RA = O*"), 
corresponding to one full season of observation. The lowest 
k bin approaches the sample variance limit as the signal to 
noise per voxel reaches ~ 1 and the array begins to image the 
largest EoR scales. The higher k bins, however, are thermal 
noise dominated at 900 hours. 

We can also follow iLidz et al] (|2008l ) and fit an ampli- 
tude and slope to In A2i(fc) in In(fc), 

In Aii(fc) = In Aii(fc = kp) -f aln(fc/fcp), (5) 

where kp is a fixed pivot wavenumber. The uncertainty on 
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the amplitude depends on the pivot wavenumber, and we 
choose kp = 0.06 Mpc"^. The uncertainly is estimated as- 
suming Gaussian statistics, and we fit directly in the 3D 
fc-space to avoid binning effects and biases. For a full sea- 
son of observation (900 hours on a primary field, 700 hours 
on a secondary), we predict a SNR of 14 on the amplitude 
and 10.9 on the slope (a) for the fiducial model. This does 
not take into account instrument downtime due to inevitable 
maintenance, nor loss of data for unforeseen reasons. With a 
more conservative observation time of 450 hours on a single 
field, we expect a SNR of 7.1 on the amplitude and 5.0 on 
the slope. Even with less than half a full observing season, 
the MWA has the potential for an EoR detection. 

This calculation does not account for systematic biases 
from calibration and foreground subtraction errors. Efforts 
are underway to un derstand these aff ects and to achieve this 
level of sensitivity l|Trott et ahlbOld ). 



3 CONCLUSIONS 

Using the proposed 128 antenna MWA, we have estimated 
the instrument sensitivity to a model EoR power spectrum, 
taking into account synthesis rotation, chromatic and asym- 
metrical baseline effects, and excluding modes that are con- 
taminated by foreground subtraction. We provide the tools 
required to calculate the MWA sensitivity for any model. 
With an optimistic full season of observation, we would ex- 
pect to detect the fiducial power spectrum amplitude with 
SNR ~ 14, and constrain the slope with SNR ~ 10.9. As 
of mid- July, construction well underway on the MWA, and 
first light is expected at the end of 2012. 
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SUPPLEMENTARY MATERIAL 

Two supplementary tables are available for electronic down- 
load. First, MWA_Antenna_Locations.txt provides the abso- 



lute coordinates of the antennas as they have been placed. 
The first column is the antenna number. The second and 
third columns are the Easting and Northing coordinates for 
a UTM projection in meters. The fourth column is the ele- 
vation in meters. All coordinates refer to the south west cor- 
ner of the antennas. The second table, obs_time_table.txt, 
provides the expected observed time per uv cell per day. 
This table is the data used to produce Figure [2] The first 
two columns are the (it, v) coordinates in wavelengths at 158 
MHz. The third column is the observation time per day for 
the corresponding uv cell in seconds. The grid size is 8.3 
m du e to the instrumental window function (|Bowman et al.l 
l2006l ). The data has been padded with zeros so the array 
can be reshaped to create a two dimensional image on a 
regularly spaced grid (see header text). 
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